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CALCULATED  PEAK  PRESSURE -DISTANCE  CURVES  FOR  PENTOLITE  AND  TNT 

ABSTRACT 

The  Kirkvood-Brinkley  shock  propagation  theory  is  used  to  obtain 
peak  pressure-distance  curves  for  Pentolite  and  TNT.  The  curve  for 
Pentolite  obtained  by  using  the  calculated  initial  pressure  and  flow 
ener©r  at  the  charge  surface  as  initial  values  is  in  excellent  agree¬ 
ment  with  experiment  for  distances  greater  than  I.75  charge  radii. 

The  curve  for  TNT  obtained  by  using  some  approximate  values  of 
the  peak  pressure  and  ener©r  at  the  charge  surface  is  in  excellent 
agreement  with  some  recent  free -air  blast  measurements.  The  curve  is 
also  compared  with  Brode's  calculated  pressure-distance  curve. 


TABLE  OF  CONTENTS 


ABSTRACT  . 

SYMBOLS  . 

INTRODUCTION  . 

THE  KIRKWOOD -BRINKLEY  THEORY . 

INi'LOKAj.’xON  OF  K-Lt^AWL^juBdvxiiKLEY  EQUATIONS* 

RESULTS  . 

ACKNW  ,’LEDGEMEWr . 

REFERENCES  .  .  ,  ,  . 

FIGURE  1  . 

FIGURE  2  . 

DISTRIBUTION  LIST  . 


PAGE 

t 

s 

6 

•7 

9 

1  7 
> 

14 

1 6 

17 

19 

20 
21 


5 


LIST  OF  SYMBOLS 


c  local  sound  velocity 

cq  sound  velocity  at  3&0°K  and  1  atm. 

E  energy 

h  enthalpy 

p  excess  pressure 

P0  ambient  pressure,  1  atm. 

r  distance 

f  f  4  *na 

w  w  t  ||»» 

u  mass  velocity 

U  shock  velocity 

p  density 

pQ  density  of  air  at  1  atm  and  300°K 
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INTRODUCTION 


The  prediction  of  damage  to  structures  caused  by  blast  waves  in 
air  requires  knowledge  of  the  peak  pressure-distance,  pres sure -time, 
reflected  pressure-time  relations  as  well  as  the  impulse-distance 
relation.  Extensive  experimental  data  exist  for  the  propagation  of 
the  blast  wave  from  bare  spheric'!  charges  detonated  in  air  and  many 
of  these  data  (for  Pentolite)  have  been  compiled  by  Goodman^-.  The 
theoretical  determination  of  the  blast  wave  paramenters  requires  the 
numerical,  integration  of  the  basic  hydro dynamical  equations  describing 

the  flow  end  until  recently  f ?y  «*•?••) ^factory  reeulte  were  available . 

2 

Brode  has  calculated  the  blast  wave  from  a  spherical  charge  of  TNT  and 

reports  the  shock  line  parameters  as  well  as  the  flow  properties  behind  the 

3 

shock.  In  a  less  extensive  calculation,  Makino  and  Shear  have  calculated 
the  flow  properties  behind  an  experimentally  "known"  shock  line  for 
Pentolite  by  the  method  of  characteristics.  The  calculation  of  Makino 
and  Shear  gives  the  flow  properties  in  the  region  bounded  by  the  shock  line, 
the  explosion  gas -air  boundary  and  a  terminal  characteristic]  hence  does 
not  include  calculations  of  impulse.  Neither  of  these  calculations 
include  shock  reflection  from  a  rigid  boundary.  Therefore,  one  must 
rely  upon  experimental  data  for  the  flow  properties  during  reflection^" 
or  estimates  of  the  shock  reflection  parameters  which  are  based  upon 
knowledge  of  the  peak  reflected  pressure  or  assumed  time  histories  of 
pressure,  (e.g.,  ref.  b). 

The  peak  pressure -distance  curve  for  Pentolite  has  been  determined 
by  extensive  experimental  data.  This  curve,  given  appropriate  initial 
data,  could  be  calculated  by  the  method  used  by  Brode  in  his  calculations 
for  TNT  and  comparisons  made  with  the  experimental  data.  The  pressure- 
distance  curve  can  be  calculated  by  the  shock  propagation  theory  of 
Kirkvood-Brinkley^  if  appropriate  initial  conditions  are  given.  The 
early  results  of  Kirkwood  and  Brinkley  showed  good  agreement  with 
experiment.  Since  that  time  the  equations  of  state  for  air  have  been 
revised,  experimental  techniques  have  been  Improved  and  great  mass  of 
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experimental  data  has  been  accumulated.  Thus,  since  the  Kirkwood- 
Brinkley  theory  does  provide  a  simple  method  of  obtaining  the  peak 
pressure-distance  curve,  a  comparison  should  be  made  vith  existing 
data.  Also,  comparisons  should  be  made  with  the  more  refined  cal¬ 
culations  of  the  type  performed  by  Brode,  when  such  calculations 
become  available. 

In  the  following  pages,  the  Kirkwood-Brinkley  theory  will  be 
described  briefly  and  used  to  obtain  pres sure -distance  curves  for 
Pentol ite  and  T*?T  under  various  initial  conditions.  The  initial  con¬ 
ditions  which  are  required  in  the  integration  are  to  be  specified  at 
(l)  the  charge  surface  or  (2)  at  a  given  fixed  distance  from  the  center 
of  the  charge. 


THE  KIRKWOOD-BRINKLEY  THEORY 
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The  shock  propagation  theory  of  Kirkwood  and.  Brinkley  is  based 
upon  a  similarity  constraint  imposed  upon  the  energy-time  curve  of  the 
shock  wave.  The  use  of  this  similarity  constraint  and  the  Hugoniot 
relations  enables  one  to  reduce  the  partial  differential  equations 
describing  the  flow  process  to  a  system  of  ordinary  differential  equa¬ 
tions  for  peak  pressure  and  shock  wave  energy  as  functions  of  distance 
from  the  explosive  charge.  Detailed  derivations  and  analysis  of  this 
propagation  theory  may  be  found  in  references  5>  6>  and  7;  hence  it 
will  be  sufficient  to  outline  the  basic  relations  involved  in  the 
approximation. 

Kirkwood  and  Brinkley  use  the  equations  of  continuity  and  motion, 
which  describe  the  flow,  in  the  form: 


where  the  variable  R  is  the  position  in  the  undisturbed  air  of  a 
volume  clement  which  has  at  tine  t  the  position  r,  thus  at  the  shock 
front  r  a  R.  Thus,  at  the  shock  front  the  equations  of  continuity 
and  motion  become 


and  the  Hugoniot  relation 


p  =  pQ  Uu 


S 
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(2) 


when  differentiated  in  the  direction  of  the  shock  path  gives  a  third 
relation,  viz., 

(&)  ♦  <!)  ■  »  fP(»  -  |(sf) 

R  t  / 1  '  R 

The  shock  front  parameters,  U;  p  and  u  are  given  as  functions  of 
p  by  the  Hugoniot  relations.  To  solve  equations  (l)  and  (3)  algebraically 
for  the  four  partial  derivations,  a  fourth  relation  among  the  partial 
derivatives  is  needed.  Kirkwood  and  Brinkely  obtain  the  deBired  fourth 
condition  by  considering  the  total  work  done  by  the  explosive  source  and 
through  their  analysis  express  the  energy-time  integral 


*oo 


r  pu  dt 


in  terms  of  the  peak  pressure-distance  curve  of  the  shock  wave  at 
distances  beyond  R,  that  is, 

*oo  oo 
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2 

r  pu  dt 
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po  R 


Ah  dR 


(4) 


where  Ah  le  the  specific  enthalpy  increment  of  a  fluid  element  traversed 
by  the  shock  wavej  the  shock  leaves  the  element  in  a  state  of  higher 
entropy  and  energy,  from  which  state  it  returns  adiabatr.cally  to  ambient 
pressure.  Hence  Ah  depends  only  on  the  peak  pressure  at  R,  i.e.,  Ah  a 
Ah  (p  ).  normalizing  the  time  Integral  in  (k)  by  the  shock  front  values 


R  PpU^  and  choosing  a  reduced  time  t  given  by 


T  = 


(5) 


where 
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in  (r2  p  u  ) 


t  a  0 
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The  energy -time  integral,  E,  becomes 


i 


l 


E(R) 


UR  v(R)n 


where 

v(R) 


and  f(R,T)  = 


f  (R,  t)  dr 
2 

r  pu 
R2  PR  “R 


(?) 

(8) 

(9) 


The  basis  of  the  Kirkvood-Brinkley  shock  propagation  theory  is 
the  approximation  of  the  integral  v  (R)  (eq.  8).  Kirkwood,  and  Brinkley 
(ref.  6)  find  as  a  satisfactory  approximation  that 


v(R) 


1 


1/5  exp 


(10) 


With  this  approximation  and  equations  (6)  and  (7)  the  desired  fourth 
condition  among  the  partial  derivatives  is  obtained.  Thus,  the  equation 


1  8  u  .  1  a  p  ,  2u 
u  ST  p  St  F 


1  S 
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2u 
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and  equations  (l)  and  (3)  supplemented  by  the  Hugoniot  relations  give 
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form 
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4(pQ/p)  +  2  (1  -  p  /p)  G 

where  ft  =  ft  (p/pq)  =  1  -  2(l  ^  - 2 - 

m(p/p0)  =  (r)  2(1  ^e)  :-qT 

0  4  ^ 


H(p/p0)  =  M 


7  +  1 


2 

°o  P 


x  *  in  (R/R  ) 
y  =  in  p/p„ 


z  n  in  Q 


,  _  B  o  _  E  jjKJ 

s  U  dp 


•  ■  >  -t¥)' 


Eoe  *1 
R_ 


=  charge  i*adlus 

p  -  loading  density  of  explosive. 

© 

and  where  7  =  1.4  is  the  specific  heat  ratio  for  air  at  500°  C  and  1 
atmosphere  pressure. 
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INTEGRATION  OF  THE  KIRKWOOD -BRINKLEY  EQUATIONS 

The  coefficients  .0 ,  M  and  N  which  occur  in  equations  (12)  are  functions 

of  the  shock  pressure  only.  These  have  been  computed  from  the  tables  of 

8 

Shear  and  Day  .  Equations  (12)  cen  then  be  solved  by  specifying  either 
p  and  E  at  the  charge  surface  or  by  specifying  p  ana  the  slope  of  the 
pressure-distance  curve  at  some  fixed  value  of  R. 

9 

Shear  has  calculated  the  initial  pressure  and  flow  energy  for 
Pentolite  detonating  in  free  air;  he  reports  the  initial  pressure  to  be 
669  atm  and  the  specific  flow  energy  to  be  1+7568  cnratm/g.  These  values 
were  used  as  initial  conditions  at  R/R^  =  1,  (the  charge  surface)  in  the 
integration  of  equations  (12)  and  the  resulting  pressure  distance  curve  is 

piTSSSiitsci  i 2,. 

Sultanoff  and  McVey10  have  performed  some  very  precise  measurements 
of  shock  distance  vs.  time  for  Pentolite.  At  a  distance  corresponding 
to  20  charge  radii  the  pressure  is  reported  to  be  10.108  atm.  and  the 
slope,  calculated  by  a  5-poiQb  Lagrangian  interpolated  formula  is 
-2.1+35-  These  values  were  used  as  initial  values  in  equations  (12)  and 
the  calculated  pressure-distance  curve  curve  is  shown  in  Figure  1. 

The  initial  conditions,  at  the  charge  surface,  for  TNT  appear  to  be 

2 

suspect.  Brode  used  an  initial  flow  energy  of  252.28  K  cal/mole  TNT 
in  his  calculations  but  because  of  inconsistencies  in  his  Figure  1  and 
Figure  16  it  it  difficult  to  estimate  the  initial  pressure.  When  the 
pressure-distance  curve  of  Brode* s  Figure  1  is  extrapolated  back  to  the 
charge  surface  the  initial  pressure  appears  to  be  between  1+80  and  490 
atmospheres,  whereas  Figure  16  of  the  came  report  shows  a  peak  pressure 
slightly  less  than  400  atm. 

Shear1*  indicates  that  he  obtained  an  approximate  value  of  515 

atmospheres  for  the  initial  peak  pressure  from  Jones  and  Hiller's  data11 

and  Brode  reports  on  initial  flow  energy  of  221.7  K  cal/mole  TlfT  that 

12 

he  calculated  from  the  work  of  G.I.  Taylor  .  Since  Taylor's  and  Jones 
and  Miller's  data  are  not  in  agreement,  consistency  is  again  lacking. 
However,  using  P  =  515  atm  and  E  =  221.7  cal/mole  TNT  at  the  charge 
surface,  equations  (12)  were  integrated  and  the  resulting  pressure- 
distance  curve  is  presented  in  Figure  2. 
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RESULTS 


The  pressure-distance  curve  obtained  from  the  Klrkvood-Brinkley 

9 

equations  (12)  using  Shear's  initial  data  shows  excellent  agreement 
with  experimental  data  except  in  the  region  of  the  charge  surface  (see 
Figure  1)  i.e.,  except  for  R/R^  between  1  and  1.75*  The  pressure- 
distance  curve  obtained  using  the  pressure  and  the  slope  of  the  experimental 
pressure-distance  curve"^  at  20  charge  radii  shows  very  good  agreement 
with  experiment  for  R/R^—  3.  At  r/r^  =  1  the  pressure  is  about  570  atm 
(see  Figure  l). 

15  11 

Figure  2  shows  some  experimental  pres sure -distance  data  * 
along  with  Brode's  curve  for  TNT  and  the  pressure-distance  curve  obtained 
from  solving  equations  (12)  with  initial  values  of  pressure  and  energy 
of  515  atm  and  221. 7  K  cal/mole  TOT,  respectively.  The  data  of  reference 
15  were  obtained  from  the  explosion  of  7.9-pound  spherical  charges  of  TOT 
detonated  in  free  air  and  the  data  rf  reference  ll  were  obtained  from 
the  detonation  of  a  20-ton  hemispherical  charge  of  TOT  exploded  on  the 
ground.  In  plotting  the  data  of  reference  14,  no  corrections  were  made 
for  shock  interaction  with  the  ground  or  the  "so-called"  reflectivity 
factors.  Brode's  curve  seems  to  describe  the  measurements  of  reference  ll 
whereas  the  calculated  Kirkwood-Brinkley  curve  reproduces  the  free-air 
measurements  of  reference  1J.  For  R/R^IO  Brode's  curve  and  the  data  of 
reference  14  can  be  shifted  to  show  agreement  with  the  Kirkwood-Brinkley 
curve  and  the  data  of  reference  15  but  this  requires  an  increase  of 
approximately  JO  percent  in  the  eners’’  used  by  Brode  in  his  calculations. 

The  difference  in  the  shape  of  the  initial  protions  of  the  two  curves 
prohibits  any  possibility  of  overall  agreement  for  the  two  curves. 

Pressure-diBtance  curves  obtained  by  the  Klrkvood-Brinkley  theory 
with  an  initial  energy  of  252.2  K  cal/mole  TOT  and  initial  pre&sureB 
in  the  range  4JO-515  atm.  did  not  differ  greatly  from  the  plotted 
curve  in  Figure  2  for  values  of  R/R^>10.  For  R/Rj>720  these  calculations 
give  pressures  higher  than  the  corresponding  pressures  of  the  plotted  curve. 

Kirkwood  and  Brinkley  calculated  a  theoretical  pressure -distance 
curve  for  cast  TOT.  (UTRC.  A-341).  They  used  p=525  atmj  and  E  =  240.6 
osl/mole  TOT  in  their  calculations j  the  resulting  curve  does  not  differ 
significantly  from  the  curve  plotted  in  Figure  2. 

lb 


The  excellent  agreement  of  the  pressure-distance  curve  for  Pentolite, 
calculated  from  the  data  of  Shear,  with  the  experimental  pressure-distance 
data  indicates  the  need  for  re -investigation  of  the  Kirkwood-Brinkley 
theory  and  further  calculations  of  the  detonation  characteristics  of 
explosives.  For,  if  initial  conditions  for  the  Kirkwood-Brinkley  theory 
can  he  prescribed  by  theory  and  the  resulting  shock  line  adequately 
describes  the  pressure -distance  relation,  this  provides  a  quick  and 
inexpensive  method  for  prediction  of  the  shock  properties  of  explosives 
without  need  of  any  experimental  data.  Thus,  the  number  of  costly  and 
time  consuming  experiments  required  to  determine  the  blast  properties  of 
newly  developed  explosives  could  be  reduced.  The  above  described  calcu¬ 
lations  could  be  performed  in  a  few  hours  and  the  results  could  be  used 
as  initial  conditions  for  a  blast  calculation  of  the  type  performed  by 
Brode.  The  Kirkwood-Brinkley  calculation  required  less  than  ten  minutes 
of  computation  time. 
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